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ABSTRACT 


Recent  thermal  decomposition  studies  of  SF^  have  led  to  interest 
in  its  principal  oxidation  product,  sulfuryl  difluoridc.  In  the  present  study 
the  thermal  stability  of  SO^  at  high  temperatures  has  been  investigated. 
Highly  dilute  SO^F^-Ar  mixtures  (  -  0.  1%)  were  shock  heated  in  a  conven¬ 
tional  1.  5"  stainless  steel  shock  tube.  The  SO^F^  concentration  was  moni¬ 
tored  as  a  function  of  time  behind  the  incident  shock  wave  by  its  infrared 
emission  at  II.  7  utilizing  a  liquid  helium  coolcdCu:Ge  detector.  The 
initial  pressure  in  the  shock  tube  was  varied  from  30  to  600  torr  and  the 
temperature  range  covered  was  1900-2300°K.  The  monitored  radiation 
was  shown  to  be  transparent  over  the  range  of  densities  employed.  Ef¬ 
fective  first-order  rate  constants  were  evaluated  from  the  lrgarithmir 
initial  slopes  of  the  radiation  decay  curves.  For  the  30-torr  data,  a  uni- 
molecular  rate  constant  fit  to  ths  data  is  k^j.  =  2.  1  x  10*  *  exp  (-39,200/T) 
sec’*.  The  data  arc  analyzed  in  the  light  of  several  modern  unimolccular 
rate  theories,  yielding  «  81  kcal  as  the  cndothcrmicity  for  the  reaction 
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I.  INTRODUCTION 


Recently  there  have  been  a  number  of  experimental  studies  of  the 

1  2 

kinetics  of  the  thermal  dissociation  of  SF^  ’  .  In  the  present  work  we 
investigated  the  stability  and  kinetics  of  dissociation  of  sulfuryl  difluoride, 
SO^F^,  which  is  (thermodynamically)  the  principal  oxidation  product  of 
SF^.  The  rate  information  was  obtained  in  the  neighborhood  of  2000°K. 

Highly  dilute  mixtures  of  50^2  in  Ar  were  shock-heated  to 
temperatures  above  which  SO^F^  is  unstable,  and  its  disappearance  as  a 
function  of  time  behind  the  incident  shock  wave  was  monitored  in  emission 
at  1 1 . 7/i . 

The  emission  intensity  data,  over  a  wide  range  of  partial  pressures 
of  SC^F^  showed  the  gas  to  be  "transparent"  in  the  bandpass  employed. 

The  time  history  data  showed  the  decomposition  kinetics  to  be  in  the  uni- 
molecular  fall  off  regime-  over  the  factor  of  20  in  pressure  covered  in  the 
experiment. 

The  kinetic  data  have  been  interpreted  employing  both  the  RRK 
theory^  and  the  recently  developed  complex  molecule  decomposition  theory 
of  Keck  and  Kalclk  ir*.  These  analyses  indicate  an  SC^F-F  bond  of  —  81  kcal/ 


mole. 


II.  EXPERIMENTAL  TECHNIQUE 


The  high  purity  1.  5  inch  stainless  steel  shock  tube  used  in  the 

5 

present  study  has  been  previously  described  .  It  ia  equipped  with  automatic 
diaphragm  changers  and  was  routinely  pumped  down  to  between  10"^  and 
10  ^  torr.  The  IR  optical  station  consisted  of  gold  coated  glass  mirrors 
which  were  used  to  focus  the  image  of  the  detector  element  (1  mm  diam. 
circle)  on  the  center  line  of  the  shock  tube  through  an  Irtran-4  window. 

A  liquid  helium  cooled  copper -doped  germanium  detector  was  used.  The 
filter,  placed  just  in  front  of  the  detector,  had  approximately  a  bell  shaped 
bandpass,  with  its  peak  transmission  at  1 1 . 75/lx  and  a  full  width  at  half 
height  of  0.43 /x.  The  whole  system  was  calibrated  with  a  400°C  blackbody 
source. 

The  50^2  was  obtained  from  Matheson  and  is  specified  99.  5% 
minimum  purity,  with  air  and  SO,  being  the  main  impurities.  The  pure 
argon  showed  no  emission  in  the  IR  bandpass  when  shocked  by  itself. 

Room  temperature  IR  absorption  spectra  of  a  1%  S02F2~99%  Ar  mixture 
were  run  between  2  and  15/i  in  a  10  cm  path  length  cell  at  pressures  of 
10  and  100  torr.  Bands  v^,  v^,  v ^  and  Vg  appeared  as  identified  by 
Perkins  and  Wilson^;  no  other  spectral  features  were  found. 

Mixtures  of  SO2F2  and  Ar  were  prepared  in  stainless  steel  mixing 
tanks  and  varied  from  0.  1%  to  1.  0%  SO2F2.  Runs  were  made  using  initial 
pressures  Pj  from  5  to  600  torr  covering  the  temperature  range  —  1500  - 
2300°K.  A  great  deal  of  this  data  was  used  to  establish  the  non -blackness 
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of  the  emission  lines  at  temperatures  between  about  1500  and  1900°K.  For 
the  kinetic  experiments,  where  knowledge  of  the  temperature  is  so  important, 
only  the  dilute  mixtures  could  be  used  so  as  to  minimize  the  temperature 
drop  due  to  dissociation.  The  kinetic  data  were  obtained  in  runs  at 
temperatures  between  1900  and  2300°K  for  the  following  4  cases:  0.1% 

^2^2  at  ^1  =  ^0,  300  and  100  torr  and  0.2%  SO2F2  at  Pj  =  30  torr; 
for  these  runs  conditions  were  essentially  isothermal. 

The  signal  to  noise  ratio  was  somewhat  of  a  problem  in  this  experi¬ 
ment;  it  was  typically  of  order  10.  The  main  sources  of  noise  were  the 
detector  and  its  load  resistor  and  the  preamplifier.  The  time  resolution 
was  set  by  the  shock  transit  time  across  the  detector  image  and  the 
electronic  rise  time  of  the  detector  circuit,  these  two  being  of  the  same 
order  of  magnitude.  The  net  effective  risetime  in  the  kinetic  experiment 
was  3.  7  and  5.  5 jx  sec  (10%  -  90%  signal  risetime),  the  faster  one  being 
used  at  the  higher  temperatures  to  allow  resolution  of  the  faster  chemistry. 
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III.  EXPERIMENTAL  RESULTS 


Two  typical  oscillograph  traces  are  shown  in  Fig.  1  representing 
(a)  the  low  temperature  (1690°K)  and  (b)  the  high  temperature  (2080°K) 
data.  At  low  temperatures  (1500  -  1900°K),  the  IR  signal  is  seen  to  be 
essentially  a  step  function,  the  rise  time  of  which  is  instrumental.  At 
the  temperature  of  the  run  shown  in  Fig.  la,  the  30^2  is  either  stable 
or  dissociates  so  slowly  that  the  IR  emission  stays  constant  following 
vibrational  excitation.  In  Fig.  lb,  the  temperature  is  sufficiently  high 
so  that  the  SO^F^  decomposes  subsequent  to  vibrational  relaxation; 
hence,  the  IR  signal  decays  as  a  function  of  time. 

These  considerations  are  illustrated  in  Fig.  2  where  the  equilibrium 

concentrations  have  been  plotted  as  a  function  of  temperature  for  the  case 

0.  1  torr  initial  partial  pressure  of  SC^F^,  e.  g.  ,  0.  1%  50^2-99.  9 %  Ar, 

Pj  =  100  torr  shock  conditions.  Two  cases  are  shown:  (1)  the  enthalpy 

of  formation  of  SC^F  was  assumed  to  be  AH^.  (298°K)  =  -  143  kcal/mole 

(S02F-F  bond  =  81  kcal)  and  (2)  AHf  (298°K)  =  -129  kcal/mole  (SC^F-F 

bond  =  95  kcal).  The  later  case  corresponds  very  nearly  to  omitting 

SO2F  completely  as  an  allowed  species  in  so  far  as  the  stability  of 

SO2F2  at  high  temperature  is  concerned.  To  keep  the  figure  relatively 

simple,  the  F-atom  concentration  was  not  plotted,  nor  were  SO,  O2,  O 

and  S  at  the  higher  temperatures  where  the  SO2  is  dissociating.  These 

equilibrium  calculations  were  carried  out  on  a  computer  using  as  input 

7 

the  thermodynamic  functions  given  by  Wilkins  for  SF^  and  its  decompo¬ 
sition  and  oxidation  p;  ducts;  however,  SO2F  was  not  considered  in 
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Fig.  1  Typical  oscillograph  traces  showing:  (a)  No  dissociation; 

Us  -  1 . 26  mm/flscc,  -  30  torr ,  SO2F2  =  1.0%,  T  =  1690°K, 
and  (b)  Dissociation;  U  =  1 . 42  mm/flsec,  Pi  =  300  torr, 

S02F2  =  0.  1%,  T  -  2080°K. 
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Fig.  2  Equilibrium  concentrations  of  SO?F2,  SC^F  and  SO2  as  a 

function  of  temperature  for  initial  shock  conditions  Pj  =  100 
torr,  SO2F2  =  0.1%.  Two  cases  are  shown,  i.  e.  ,  SO2F 
heats  of  formation  of  -143  and  -129  kcal/mole,  corresponding 
to  SO2F-F  bond  strengths  of  81  and  95  kcal,  respectively. 


Ref.  7.  As  shall  be  seen  below,  the  present  experiments  indicate  an 
SO^F-F  bond  strength  of  ~  81  kcal,  and  hence  at  intermediate  temperatures 
it  becomes  a  major  specie.  The  thermodynamic  properties  of  the  SO^F 

g 

were  estimated  for  us  by  Wilkins  . 

It  would  have  been  highly  desirable  to  measure  the  SO^F-F  bond 
strength  directly  in  the  present  work  by  monitoring  the  equilibrium 
concentration  of  SO^F^  in  the  sensitive  temperature  range  between  1600 
and  2000°K.  However,  shock  attenuation  and  boundary  layer  growth 
effects  make  the  temperature  uncertain  by  up  to  several  hundred  degrees 
by  the  time  equilibrium  is  established  (~  100  jisec  in  the  present  experi¬ 
ment).  Inspection  of  Fig.  2  shows  that  this  magnitude  of  temperature 
uncertainty  doe.;  not  allow  determination  of  the  bond  strength  to  the  desired 
accuracy  by  this  method. 

With  the  present  apparatus,  meaningful  kinetic  data  could  not  be 
obtained  above  about  2300°K  due  to  insufficient  resolution.  Below  about 
1900°K,  both  the  small  amount  of  decomposition  and  the  slow  rate  of 
decomposition  prohibited  kinetic  measurements  as  the  effects  of  shock 
attenuation  were  noted  to  become  dominant  under  these  conditions. 

The  data  obtained  in  the  temperature  range  1500  -  1900°K  were 
used  to  demonstrate  the  transparency  of  the  emitting  gas.  For  a  thin  gas, 
the  intensity  in  a  band  is  given  by 


I  =  a  /  A^  jsO^F^  watts/cm^-ster  , 


(1) 


where  a  is  the  integrated  band  intensity  in  cm"^  atm  is  the1  value 

of  the  blackbody  function  at  the  radiating  gas  temperature  in  watts/cm^- 
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ster-cm,  I  is  the  optical  path  length  in  cm,  A  is  the  wavelength  in  cm 


and  SO^F^  is  the  concentration  of  SO^F^  in  amagats; 


KM 


f  . 


12) 


where  P{  is  the  initial  shock  tube  pressure  in  torr,  P^/Pj  1*  (be  density 
ratio  across  the  shock,  Tj  is  the  initial  temperature  of  the  test  gas  and 
f  is  the  fraction  of  SO^F^  in  the  test  gas  mixture.  Equations  (1)  and  (2) 
indicate  that  for  a  thin  gas 


n~  *  f  (3) 

*1 

The  plateau  values  ol  the  IR  signals,  normalized  by  the  pressure, 
have  been  plotted  against  the  percent  of  S02^2  in  Fig.  3.  Each  one  of 
these  data  points  Is  an  average  of  from  3-10  separate  runs  covering  a 
temperature  range  of  1500  -  I900°K.  Within  the  scatter  in  the  data,  the 
temperature  dependence  over  this  range  was  negligible.  The  data  fall  on 
a  line  of  unit  slope  on  the  log -log  plot  indicating  agreement  with  Eq.  (3) 
and,  hence,  that  the  emission  lines  are  not  black. 

It  was  important  to  show  that  under  the  experimental  conditions, 

Eq.  *3)  is  followed.  This  unambiguously  shows  that  the  radiation  signal 
is  proportional  to  the  SO^F^  concentration,  thus  allowing  interpretation  of 
the  kinetic  data.  Before  leaving  the  data  shown  In  Fig,  3,  however,  it  is 
of  some  value  to  use  it  to  obts'n  an  estimate  ofihe  integrated  band  intensity 
a  .  We  note  that 
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FRACTION  SGgf^ 


Fig.  J  Peak  ll.7fi  radiation  signals,  normalised  with  the  initial 
shock  tube  pressure  plotted  vs.  parcant  of  SO»F2  in  the 
ft  a  •  mixture.  Each  data  point  is  the  Average  of  many  runs. 
1M>0£T  <  I900«K. 
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t*4*K 


a 


(4) 


jJ-  760/.V  . 

P,  k  p,  T, 


*nd  that  l/Pj  if  obtainable  from  Fig.  3  by  u«e  of  the  btackbody  calibration 
factor  J  and  the  bandpass  of  the  filter  AX  *  4.  j  x  10"^  cm: 


i 

*7 


SR  J  AX 


N 


and  J  a  ■ 


X<B) 

57“ 


(5) 


where  SR  and  SR  are  th»  detector  signals  for  a  run  and  calibration,  respec¬ 
tively,  and  N^ipj  it  the  value  of  the  biackbody  function  at  the  calibration 
source  temperature.  From  Eqs.  (4)  and  (3)  and  the  data  of  Fig.  3  we  find 
o  *  500  cm*^  atm"*. 

It  should  be  pointed  out  that  the  above  value  of  o  is,  at  best,  a* 
approximate  one.  At  room  temperature^  the  (A ^ )  band  at  1 1. 84 
(848  cm"*)  is  just  about  encompassed  by  the  bandpass  filter  used  in  the 
present  experiment,  and  the  i»g  (B^l  band  at  II.  34  (885  cir."*}  is  just 
about  excluded.  However,  at  elevated  temperatures,  both  bands  must 
spread  in  wavelength  so  that  the  experimental  bandpass  must  exclude  some 
of  the  ^  and  include  some  of  the  ^  radiation. 

The  high  temperature  data  were  used  to  obtain  the  dissociation  raftc 
constant.  For  the  reaction 


SO^F^  *  products  Rl 

an  experimental  (effective)  unlmolecuiar  rate  constant  k^  is  defined  by: 


[SOjFj]  <J. 


9  k 


eff 


(4) 
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where  (SO-jF^J  is  the  concentration  of  SO^F,  ir  particlcs/cc  and  t  is  the 
particle  time  in  .  Noting  that  S  *  [SO^F  ^J  ,  where  S  is  the  detector 

f*2 

output  signal,  and  that  t  ■  - —  t.  where  t.  is  time  in  the  laboratory 

P  P|  u 

coordinate  system,  we  have 


1  dS  pl  1  dS  . 

S  3t^“  ’  *  5  3t^  ’  eff  ‘ 


<7) 


For  the  kinetic  runs  ^  ^  independent  of  mixture  or  shock  speed. 

Such  a  logarithmic  derivative  can  be  taken  any  place  in  the  time  history 
provided  that  the  decay  is  truly  exponential.  Non«exponential  decay  is 
observed  at  long  times  when  attenuation  effects  become  non-ncgligiblc, 
and  at  low  temperatures  where  the  backward  reaction  becomes  important. 
dS/dtj^  is  indicated  on  the  oscilloscope  trace  shown  in  Fig.  lb. 

The  kinetic  data  are  presented  in  Fig.  4  as  an  Arrhenius  plot  of 
k^  vs.  reciprocal  temperature.  The  four  straight  lines  are  least  square 
fits  to  the  four  sets  of  data  assuming  an  Arrhenius  type  temperature 
dependence.  For  the  30  torr  case 


kcf(  •  2.  I  x  1011  exp  (-  39,200/T)  sec'1  ,  (8) 


which  corresponds  to  an  activation  energy  of  E  =  78  kcal/mole. 

The  density  corresponding  to  the  data  for  the  extreme  P|  cases  are 
indicated  in  the  figure  as  P^/Pq>  where  pQ  is  standard  density,  ltisofsome 
interest  to  compare  the  rate  of  decomposition  of  SO2F2  with  its  pa  rent  com  pound, 
SF^.  This  Is  done  in  Fig.  4at2000°K  where  the  range  of  ke^  for  SF^  decompo¬ 
sition  is  shown  for  the  same  density  range  as  encompassed  by  the  present  data. 
This  information  was  obtained  by  cross  plotting  the  data  of  Dott  and 
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Fig.  4  Effective  rate  constant  for  SO^F^  decomposition  by  Av* 

determined  from  the  initial  slopes  plotted  vs.  reciprocal 
temperature.  The  straight  lines  are  least  square  fits  to 
the  data  taken  at  the  four  initial  pressures  of  P|  a  30,  100 
300  and  600  torr.  Shown  for  comparison  is  the  rate  constant 
for  SF^  decomposition  over  the  same  density  range. 
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Jacob*.  *  It  is  seen  that  SF^  decomposes  about  »0  times  faster  than  does 
its  oxidation  produev,  SO^F^. 

The  limits  of  Pj  and  T  were  dictated  by  the  experiment;  low  partial 
pressures  of  SO^F^  were  required  to  maintain  isothermal  conditions,  at 
low  temperatures  the  slow  rate  and  stability  of  50^2  led  to  flat  signals, 
at  high  temperatures  rise  time  limitations  prevailed;  at  low  pressures  the 
data  were  noise  limited,  and  at  high  pressures  severe  shock  attenuation 
became  dominant. 
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IV.  COMPARISON  WITH  THEORY 


t 


The  kinetic  data  have  been  analysed  within  the  framework  of  two 
complex  molecule  dissociation  theories,  i.e.,  those  of  Rice-.  Ramspcrger 
and  Kassel  (RRK)*  and  Keck  and  Kalelkar*.  The  well  known  RRK  theory 
has  in  it  a  number  of  adjustable  parameters.  Figure  5  shows  the 
sensitivity  to  these  parameters  when  the  experimental  data  is  compared 
to  the  theory.  In  this  figure  kd/k—  is  plotted  against  logjQ  (A/A-4, 
where  krf  is  the  theoretical  dissociation  rate  constant  which  Is  normalised 
by  k—  -  the  limiting  high  pressure  rate  constant.  A  is  the  frequency  factor, 
and  A  is  the  collision  efficiency.  The  collision  frequency  u  is  given  by 


•  ’°2  <7jrrr 


j>/2  p 


(9) 


where  a  is  the  mean  collision  diameter.  |i  is  the  reduerd  mass  for  the 
collision  pair,  and  P  is  the  total  pressure  of  the  system. 

The  hi^h  pressure  limit  to  the  rate  constant  is  given  by 

k^  A  exp  ( -  E/kT)  .  (10) 

where  E  is  the  energy  of  the  bond  which  is  broken  when  the  reactant 
molecule  dissociates. 

The  rate  constant  is  given  by 

kd/k^  1  (S,  E/kT,  log  A/Aw)  .  (11) 

where  1  is  the  Kassel  integral,  and  S  is  the  number  of  classical  harmonic 
oscillators  that  couple  to  the  dissociating  bond,  this  usually  being  about 
1  Z  of  the  total  number  of  vibrational  modes.  Tables  of  this  function 
have  been  prepared  by  G.  Emanuel. 
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Fig.  5  Sensitivity  of  RRK  thoorcfic.il  rate  constant  to  the  theory 
parameters  and  comparison  with  data  at  2000  and  2200°K. 
The  parameters  S  and  A  were  varied  for  assumed  bond 
strengths  of  hi  and  95  keal. 
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Hence,  the  parameters  which  are  needed  are  a,  A  ,  A,  S  and  E. 
The  molecular  diameter  of  SC^F^  is  not  known;  a  value  of  n  =  4  R.  for 
the  Ar  -  SO^F^  collision  pair  was  chosen  and  cannot  be  in  error  by  a 
significant  amount.  Following  the  work  of  Bott  &  Jacobs^,  A  was  set 
equal  to  l/4.  A,  S  and  E  were  adjusted  to  obtain  the  best  fits  with  the 
data.  It  was  assumed  that  the  reaction  mechanism  was 

SOzF2  f  Ar  -  S02F  +  F  +  Ar  ,  R2 

as  this  clearly  has  the  lowest  activation  energy  of  the  several  conceivable 
SC>2F2  decomposition  paths.  Two  values  of  E  were  considered,  81  and 
95  kcal.  The  latter  was  chosen,  rather  arbitrarily,  based  on  the  average 
S-F  bond  energy  of  86  kcal  in  S02F2.  The  value  of  81  kcal  was  estimated 
by  Wilkins®. 

As  is  seen  in  Fig.  5,  F  =  95  requires  S  =  7  or  larger  at  T  = 

2000°K;  the  best,  but  not  very  satisfactory,  fit  being  made  with 
14  -1 

A  =  3  x  10  sec  .  The  E  =  81  case  can  be  only  poorly  fit  as  S  =  5  with 

12  13 

A  =  4  x  10  but  is  nicely  fit  at  S  =  4  with  A  «  1  x  10  .  Similarly,  a 

comparison  of  theory  and  data  at  2200°K  shows  that  E  =  81,  S  =  4  and 

1 3 

A  =  1  x  10  is  a  good  fit  (see  Fig.  5). 

In  Fig.  6,  all  the  data  is  compared  to  the  RRK  theory  in  the  k^, 

1/T  plane  for  the  best  results  for  the  cases  E  =  95  and  81  as  shown  in 

1  3 

Fig.  5.  Here  it  is  seen  that  the  case  E  =  81,  S  =  4,  A  =  l  x  10  fits 

the  data  adequately,  while  the  slope  of  the  E  =  95  case  is  too  steep.  A 

13  -1 

frequency  factor  of  1  x  1  0  sec  is  quite  typical  for  reactions  of  the 
type  studied  in  the  present  work  and  S  =  4  is  "reasonable"  for  the  S02F2> 
which  has  9  vibrational  degrees  of  freedom. 
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The  theory  of  Keck  &  Kalelkar^is  similar  to  that  of  RRKN1*  in 
its  basic  assumptions  but  it  differs  in  detail  and  the  results  are  in  a 


simple  form  v-'hich  facilitates  computations.  A  successful  correlation 

of  for  21  molecules  having  3-6  atoms  has  been  obtained  with  this 
4 

theory  .  The  theory  has  three  adjustable  parameters,  E,  ^  and 
E  has  the  same  interpretation  here  as  in  the  RRK  theory,  is  the 

effective  energy  transfer  cross  section  and  o '  ^  is  the  effective  spontaneous 
decay  cross  section.  These  cross  sections,  within  the  framework  of  the 
theory,  are  not  completely  adjustable  but  *  3  and,  for  the  present 
reaction  where  one  bond  is  broken  and  an  atom  is  one  of  the  products, 


A  comparison  between  the  data  and  the  Keck  theory  at  2000°K  is 

/  P2  °1A 

shown  in  Fig.  7,  where  log  k^/o^  i*  plotted  against  log- —  - - .  It  is 

seen  that  for  the  case  o^  =  3  R,  =  1  R,  as  called  for  by  the  theory, 
the  comparison  between  theory  and  experiment  is  very  poor,  and  E  >  95 


kcal  is  required.  Better  fits  can  be  made  at  E  =81,  but  O^  =  o^  =  0.  04 
for  this  case,  which  is  not  allowed  by  the  model.  In  the  k^  vs.  l/T 


plane,  the  Keck  theory  does  not  fit  the  data  well  either. 
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Fig.  7  Sensitivity  of  Keck  theoretical  rate  con  tant  to  the  theory 
parameter*  and  comparison  with  data  at  2000'^.  The 
parameters  >i|^  and  o^j  were  varied  for  assumed  bond 
strengths  of  HI  and  95  keal. 
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